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Abstract 

This paper studies foot bone geometrical shape and its mass distribution and establishes an assessment 
method of bone strength. Using spiral CT scanning, with an accuracy of sub-millimeter, we analyze 
the data of 384 pieces of foot bones in vivo and investigate the relationship between the bone's external 
shape and internal structure. This analysis is explored on the bases of the bone's center of mass and its 
centroid of shape. We observe the phenomenon of superposition of center of mass and centroid of shape 
fairly precisely, indicating a possible appearance of biomechanical organism. We investigate two aspects 
of the geometrical shape, (i) distance between compact bone's centroid of shape and that of the bone 
and (ii) the mean radius of the same density bone issue relative to the bone's centroid of shape. These 
quantities are used to interpret the influence of different physical exercises imposed on bone strength, 
thereby contributing to an alternate assessment technique to bone strength. 

Introduction 

The structure and function of bone largely depends on its mechanical and biological environments [1]- 
[5]. When bone is bearing external force, bone tissue undergoes adaptive changes such as physiological 
activities of remodeling and reconstruction [6]. These changes emerge as those of external shape and 
internal structure so as to maximize its potential to bear external load [7]-[9]. Bone's adaptive changes 
to the external force are regarded as an optimization process [Ulinilll] which observes some special law. 
Current research results, however, did not provide us with convincing explanations to this special law. 

Bone strength is an important index to diagnose osteoporosis |12j . That justifies why bone strength 
has always been an issue in biomechanical research. Bone strength is determined by factors such as 
bone external shape and internal structure. How to measure bone shape and structure accurately largely 
depends on equipment. As a result, the advanced measurement equipment affects bone biomechanical 
research. With the development of CT scanning technology, more reliable modeling methods emerged. 
Many mechanical models, both macroscopically and microcosmically, depict the relationship between 
bone density, geometric shape, internal structure and bone strength |13H16| . This has made a quantitative 
assessment of bone strength more reliable. But in clinical medicine, no other quantitative assessment 
index has been more widely used than bone mineral density, which actually cannot objectively reflect the 
real condition of bone strength |14[|16lfT7] . It is imperative to explore a quantitative and practical bone 
strength assessment index. 

External force can change bone shape and structure |18| and in turn, bone shape and structure 
affect the result of external force. When we do research on the interaction of one body with another, 
bone's centroid of shape (COS) and center of mass (COM) have been considered as two basic physical 
quantities. Concerning the above-mentioned two issues, this paper focuses on selecting COS (shape) and 
COM (structure) to study the variation laws of bone in its adaptation process to external force and based 
upon these two quantities, we propose a couple of bone strength assessment indexes. 
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Materials and methods 

Previous studies P"5^ET| have shown that the isotropic materials' positions of COM and COS coincide 
since these materials form a closed three dimensional region. On the other hand, there exists enough 
experimental evidence that signatures the bone to be anisotropic in which case the COM and 

COS will superpose only under specific conditions. The broad understanding is still limited and suffers 
the lack of comprehensive theoretical model that shows a qualitative relationship between the position 
of COM and COS in anisotropic materials. The anisotropy of bone acoustic properties may also have a 
negative impact on imprecision. For example, the orientation of the calcaneus, not only around the leg 
axis but also around the heel-toe axis, must be controlled. Small movements of the leg or of the knee 
may slightly rotate the calcaneus around its long axis, resulting in undesired variations of the observables 
and positioning errors. Whereas several devices offer means of minimizing these variations, wc choose CT 
scanning of foot bone to assess bone strength in vivo and examine the relationship between the positions 
of COM and COS in anisotropic materials. 

The test equipment is Brilliance 64-slice Scanner by Philips, Netherlands, provided by Image Process- 
ing Center of Zhujiang Hospital. Scan settings are: frame bone tissue; power: 120kv; pixel size: 0.50mm; 
layer distance: 0.50mm. The scanning is conducted along both feet transect, from top to bottom. 

Altogether, we have collected data of 384 pieces of bone - both from the voUeyballers (with aver- 
age height, weight and age of 183.94 ± 3.90cm, 69.80 ± 5.20kg and 21.88 ± 0.99yrs, respectively) and 
wrestlers (with average height, weight and age of 168.00 ± 5.68cm, 65.52 ± 5.16% and 21.00 ± 2.78yrs, 
respectively), i.e. 32 pieces of 12 types of bones: calcaneus(Ca.), talus(Ta.), navicular(Na.), cuboid(Cu.), 
lateral cunciform(La.), intermediate cuneiform(In.), medial cunciform(Mc.), first mctatarsal(Fi.), second 
metatarsal(Se.), third metatarsal(Th.), fourth metatarsal(Fo.) and the fifth metatarsal(Fif.). The num- 
ber of all of the differential elements of volume is 71,836,054. The subjects are male volleyball players 
from our institute and male wrestlers from Provincial Sports School. It has been confirmed before the 
test that every subject has been trained as a professional player for more than five years. Before the test, 
each subject's medical history was inquired and all the subjects were x-rayed to exclude subjects with 
diseases such as foot pathological change, deformity or injury to make sure that their physical conditions 
meet the requirements of the test. This study has been carried out according to the existing rules and 
regulations of our institute's Ethnic Committee. 



Results 

Following [531[5S], we separate foot bone to calculate the volume, surface area and average bone density 
(Table [1]). Due to their significant differences in height and weight, there exist statistically significant 
differences between the voUeyballers' and wrestlers' foot bone volume and surface area. However, the 
acquired differential element of volume (dv) is a constant. To calculate bone's COM and COS, we consider 
n number of infinitesimal volume elements dv with density p. Choosing the orthogonal coordinate system 
to locate the position of elementary volume dv, we use 
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to evaluate the COM and COS of the bone, respectively. Setting the COS as the origin, the position 
of the bone's COM relative to the COS can be calculated by calculating the spatial separation between 
any two points in three dimensions (Fig. [T]). We notice that about 99.22%; 99.22% and 99.74% of the 
distance between the foot bone COM and COS is within the range of ztlmm from X; Y and Z axes. 
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Table 1: Foot bone volume, surface area and bone density (MeaniSD) 





Wrestler 

Volume(cm"') 


Area(cm^) 


Density((;/cm'^) 


Volleyballer 

Volume(cm"^) 


Area(cm^) 


Density((7/cm'') 


Ca. 


71.01 ±8.46 


107.39 ± 8.83 


1.47 ±0.04 


83.94 ± 6.05 


120.70 ± 5.56 


1.49 ± 0.05 


Ta. 


38.30 ±4.33 


71.38 ±5.41 


1.63 ±0.04 


43.87 ± 3.33 


80.11 ±5.97 


1.65 ± 0.04 


Na. 


11.45 ± 1.39 


31.21 ±2.73 


1.56 ±0.04 


13.44 ± 1.51 


34.78 ± 3.00 


1.58 ±0.05 


Cu. 


13.87 ± 1.61 


33.14 ±2.77 


1.46 ±0.04 


15.09 ±2.69 


35.24 ± 4.78 


1.47 ±0.05 


La. 


5.91 ±0.69 


19.09 ±1.56 


1.51 ±0.04 


6.79 ±0.61 


20.99 ± 1.28 


1.53 ±0.06 


In. 


4.43 ±0.66 


15.69 ±1.56 


1.59 ±0.04 


5.20 ±0.44 


17.56 ± 1.00 


1.64 ±0.06 


Me. 


10.76 ± 1.48 


28.60 ± 2.73 


1.52 ±0.03 


12.20 ± 1.04 


31.02 ± 1.91 


1.58 ±0.05 


Fi. 


16.94 ±2.23 


44.90 ± 3.89 


1.62 ±0.05 


20.93 ±2.25 


51.94 ±3.47 


1.65 ±0.05 


Se. 


9.01 ± 1.29 


33.72 ±3.25 


1.73 ±0.07 


11.65 ±0.77 


40.56 ± 2.08 


1.76 ±0.08 


Th. 


7.72 ±0.58 


30.23 ±1.60 


1.70 ±0.05 


8.99 ±1.07 


34.50 ±2.61 


1.68 ±0.07 


Fo. 


7.47 ±0.78 


28.80 ±2.19 


1.66 ±0.04 


8.88 ±0.92 


32.97 ±2.05 


1.66 ±0.05 


Fif. 


8.83 ± 1.09 


30.92 ± 2.64 


1.72 ±0.05 


9.65 ±1.07 


33.73 ± 2.49 


1.71 ±0.05 



respectively. With the pixel size of spiral CT scanning at O.SOmm and its layer distance at O.SOmm, the 
measured spatial parameters suggest a superposition of the foot bone's COM and COS, which is collected 
and displayed in Fig. [TJ 




Figure 1: Relationship between the positions of foot bone's COM and COS. 



To explore this phenomenon of superposition in tissues, we separate the bone in terms of marrow, 
spongy bone and compact bone by applying the CT value of the corresponding bone tissue relationship 
Following the above procedure, we calculate the COM and COS of three tissues of each bone. 
The corresponding relationship is shown in Fig. [5] A similar behaviour of the distance between tissues' 
COM and COS is observed indicating the presence of superposition of the positions of COM and COS of 
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Figure 2: Relationship between COM and COS of tissues of marrow, spongy bone and compact bone (W 
refers to wrestlers and V volleyballers.) 



tissues of marrow, spongy bone and compact bone. Whether the superposition observed from the above 
analysis is a genuine destined phenomenon or an accidental event can be resolved from what follows. 

Discussion 

Having explored the superposition, the results can be used to extract the biomechanical significance of 
this principle. To this end we conduct an analysis of the moment of inertia, a physical quantity which is 
related both to the bone's geometric shape and its mass distribution. Using the optimal function 

min ipipc) = ^ TOi - Xb)^ + {vi - ybf + {zi - Zbf^ , (1) 

where Pb{xb, Vb, Zb) and Pi{xi, yi, Zi) refer to the positions of bone's COM and any position of dv relative 
to CT image coordinate, we differentiate the right side of Ec^. [T]to obtain 

^ 2m., {\xi - Xb\ + \yt - yb\ + \zt - Zb\) . 

For a closed and continuous three-dimensional geometric body, the series 

converge. Thus in the limit mi — 0, using the Abelian Theorem [31], the above equation becomes 
^TO,; {\xi - Xb\ + \yt - yb\ + \zi - Zb\) = ^rn^^ (|a;j - Xf,| + \y, - yb\ + \zi - Zb\) . 

Suppose 

rrii ^ {\xi - Xb\ + \yi - yb\ + \zi - Zb\) = 0, 
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then according to Parallel Axis Theorem, the point (xb, yb, Zb) must be the COM, and 



^ \xi - Xb\ = 0, ^ \yi 



yb\ 0' X! ~ ^• 



Since 



= nxb, y^yfe 




Zb = nzb, 



therefore 



xb = xc, yb = yc, Zb = Zc. 



This implies that 



Pc =Pb 



(2) 



Eq. [2] suggests that the optimal program for the heterogeneous bone tissue mass distribution should 
be the superposition of its COM and COS. This coincidence can avoid the emergence of eccentric force, 
enhance the bone's potential to bear the external force, optimize the bone's internal structure and then 
achieve the goal to use the minimal material to maximize its function f22j . Consequently, in order to get 
adapted to the mechanical environment, the bone undergoes continuous functional adaptive changes [H 
llOllllj. In this process, the bone always observes the optimized principle of superposition between the 
bone's COM and its COS. 

Experimental results show that the strength of compact bone of bone in vitro is many times more 
than that of the spongy bone [ISldn]- That means that the bone strength is mainly determined by 
compact bone. Fig. [2] indicates that the COM and COS of compact bone are in accordance with Eq. [21 
The position of compact bone's COS relative to the bone's COS will affect the external force. However, 
will the external force change the relationship of their positions? We calculate the position of compact 
bone's COS relative to the bone's COS. Set the coordinate origin at the bone's COS, the relationship of 
positions of the compact bone's COS relative to the bone's COS could be gained. See Fig. [3l 

When the compact bone's COS moves towards the bone's COS, what happens to the mass distribution 
of bone tissue? In order to explore this issue, we propose a concept of mean distribution radius (MDR) 
of same density tissue relative to its bone's COS and explore the relationship between the MDR and the 
density. The relationship between twelve types of foot bone density and the MDR relative to the bone's 
COS is shown in Fig. IH 

When a voUeyballer takes off to spike, the braking movement has a great impact on the calcaneus. In 
Table[Tl the volleyballers' calcaneus density is larger than that of the wrestlers. But there is no statistically 
significant difference {p = 0.1887). In Fig.[3l the distance of the volleyballers' calcaneus compact bone's 
COS to the bone's COS is shorter than that of the wrestlers and it has a statistically significant difference 
(p = 0.0178). In Fig. HI the MDR of the volleyballers' calcaneous relative to the bone's COS begins to 
be larger than that of the wrestlers from the density of compact bone on; especially when comparing this 
with the results from the marrow and spongy bone tissues, this difference is outstanding. The wrestlers' 
fierce body combats carry great strength to their fifth metatarsal from the front, rear, left and right. 
Table [T] suggests the wrestlers fifth metatarsal density is greater than that of the volleyballers, but again, 
there is no statistically significant difference {p = 0.4231). In Fig. [3l the distance of the wrestlers' fifth 
metatarsal compact bone COS to the bone's COS is shorter than that of the volleyballers and there is a 
big difference (p = 0.0098). In Fig. 21 the MDR of the wrestlers' fifth metatarsal relative to the bone's 
COS begins to become bigger from the density of compact bone on than that of the volleyballers. The 
difference grows especially when approaching their peak values. In both sport events, the impact strength 
to the first metatarsal is obvious. From Table [U we can see that in the similar morphological first to 
fifth metatarsal, the lowest density goes to the first metatarsal, which does not sound very reasonable. 
(Factors such as volume and joint segmental area might have affected bone density, which might as well 
mean the limitation of bone density assessment index.) In Fig. [Sj the distance of both athletic groups' 
first metatarsal compact bone COS to the bone's COS is the shortest. In Fig. |4j comparing with the 
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Figurc 3: Relationship between the compact bone's COS relative to the bone's COS. Calculate each 
individual bone's compact bone COS and the bone's COS. Calculate the distance between the compact 

bone's COS and the bone's COS. Apply r = J2i (^V i^i ~ ^cY + ~ Vci^ + (^i ^ -^c)^^ to calculate 
all bone tissues' mean distribution radii relative to the bone's COS and use dcompact / ^ to standardize the 
distance between the compact bone's COS and the bone's COS. 



second to the fourth metatarsal, the MDR of both athletic groups' first metatarsal relative to the bone's 
COS begins to enlarge. 

The above results and analyses suggest that greater strength sports activities can enlarge the bone 
density peak values (30l[31]. The idea that the impact strength can add bone mass significantly [321 - 134] 
has been verified by the voUeyballers' and wrestlers' foot bones. This is the result of bone's adaptation 
towards mechanical environment. When assessing the above results by the two indices, i.e. the distance 
between the compact bone's COS and the bone's COS and the MDR relative to the bone's COS, the 
results are consistent with the foot bone feature of load from both athletic groups. 



Conclusion 

This research studies the relationship between the positions of the COM and COS of foot bone in vivo by 
an analysis of sub-millimeter accuracy. The phenomenon of superposition of foot bone's COM and COS 
has been discovered. This phenomenon does not exist only in the relationship between the positions of 
bone's COM and COS but also in the relationship between the COM and COS of tissues such as marrow, 
spongy bone and compact bone. We establish the relationship between the tissues' masses and their 
positions, with which we analyze and illustrate that the phenomenon of superposition of bone's COM 
and COS is caused by the approximate symmetry of distribution of the bone tissue relative to the bone's 
COS. We define an active function and determine an optimal condition and prove the superposition of 
foot bone's COM and COS to be an optimized distribution program. 

We propose the concept of distance between the compact bone's COS and its bone's COS and discover 
the relationship between the distance (of the compact bone's COS and its bone's COS) and the loading 
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Figure 4: Relationship between the bone's density and the MDR relative to the bone's COS. We calculate 
the MDR to the bone's COS by applying fj — ^ J^^' (^\/ i^i ~ ^c)"^ + iVi — Vc)^ + (-^i ^ ^c^^ , calculate 

the MDR of the bone tissue relative to bone's COS by f = ^ J2i (^V i^i ~ ^<^Y + {Vi ~ VcY + {zi ~ -Zc)^^ 
and then standardize the MDR relative to the bone's COS by rj/f. 

type. This relationship is represented by the phenomenon that the impact strength has made the compact 
bone's COS move towards the bone's COS. This movement symbolizes a functional adaptation of bone 
in its structure. The physiological activity of the middle aged and seniors is mostly reconstruction [35] . 
When their bone masses are gradually decreasing, it is essential to look into the possibility of whether 
physical exercises can diminish the bone loss and change the movement's direction. This is meaningful 
and worthwhile research. 

We have studied the relationship between the bone's density and the MDR of foot bone relative to the 
bone's COS and we have observed the obvious difference that various exercises exert on compact bone's 
MDR. Since MDR has much to do with bone strength, the MDR relative to the bone's COS and the 
distance between compact bone's COS and the bone's COS can assess the bone strength more reliably. 

In conclusion, in the physiological activities of foot bones' adaptation to their external mechanical 
environment, they always observe the optimized distribution law of superposition between the COM and 
COS. By changing the distance between the compact bone's COS and the bone's COS and by changing 
bone tissue MDR to change bone strength, the bone can adapt to meet the requirements of muscle 
strength and training load to ensure the normal activities. The movement of compact bone's COS and 
the change of MDR all signal the adaptation of bone. 
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